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Upcoming Sunyaev-Zel'dovich surveys are expected to return ~ 10 4 intermediate mass clusters 
at high redshift. Their average masses must be known to same accuracy as desired for the dark 
energy properties. Internal to the surveys, the CMB potentially provides a source for lensing mass 
measurements whose distance is precisely known and behind all clusters. We develop statistical 
mass estimators from 6 quadratic combinations of CMB temperature and polarization fields that 
can simultaneously recover large-scale structure and cluster mass profiles. The performance of these 
estimators on idealized NFW clusters suggests that surveys with a ~ V beam and IO^iK' noise in 
uncontaminated temperature maps can make a ~ IOct detection, or equivalently a ~ 10% mass 
measurement for each 10 3 set of clusters. With internal or external acoustic scale i?-polarization 
measurements, the ET cross correlation estimator can provide a stringent test for contaminants 
on a first detection at ~ 1/3 the significance. For surveys that reach below 3/LtK', the EB cross 
correlation estimator should provide the most precise measurements and potentially the strongest 
control over contaminants. 

PACS numbers: 



I. INTRODUCTION 

Upcoming surveys for clusters utilizing the Sunyaev- 
Zel'dovich (SZ) effect in the Cosmic Microwave Back- 
ground (CMB) as a detection technique hold the promise 
to measure the properties of the dark energy to high pre- 
cision with ~ 10 4 clusters at redshifts out to z ~ 1. The 
mean mass and variance of the sample needs to be known 
to an accuracy comparable to that desired for the dark 
energy equation of state to not be a limiting factor. Since 
the SZ effect is sensitive to the temperature weighted 
baryon content of the cluster, it does not directly probe 
the total cluster mass in a model independent way. 

Fortunately, the same survey which identifies the clus- 
ters in the SZ effect can potentially also constrain their 
masses through gravitational lensing of the CMB. Uti- 
lizing the CMB as a source is also appealing in that its 
distance is both well determined and sufficiently far to 
probe even the highest redshift clusters in the sample. 

Studies of gravitational lensing of the CMB by clus- 
ters 

EL S S, S S S, typically focus on making de- 
tailed measurements of individual high mass clusters or 
on statistical reconstructions of average cluster proper- 
ties. The CMB suffers in the former case when compared 
with galaxy weak lensing because the background image, 
in this case the primary CMB, is a Gaussian random field 
whose properties must be separated from the cluster and 
its emission. 

This disadvantage will be substantially reduced in the 
latter case where the goal is to measure statistical proper- 
ties rather than unique objects, provided that large sam- 
ples of clusters are available for the analysis. Surveys 
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such as SPT and ACT should provide samples of ~ 10 4 
intermediate mass clusters at high redshift, so that a pre- 
cise and unbiased statistical measurement may become a 
realistic possibility in the near future. The CMB even 
possesses some advantages over galaxy lensing, since the 
source statistical properties and redshift are extremely 
well determined, and because lensing of the tempera- 
ture and polarization provide strong consistency checks 
against possible contamination by cluster emission. 

Techniques which use CMB lensing to reconstruct the 
lensing con verg ence, k, have been available for some 
time [1, ID, llOj ] and although these were initially in- 
tended as probes of large scale structure, they can in 
principle recover the projected mass regardless of its 
source. This certainly includes reconstructions of large 
galaxy clusters and associated statistical quantities, such 
as the cluster-convergence correlation function. On clus- 
ter length scales, this is essentially an average density 
profile and mass measurement (e.g. [TJ E3)- We will 
focus on estimators that are quadratic in the observed 
CMB fields [1, Q, since these can be implemented us- 
ing fast algorithms, and can eliminate contamination by 
isolating pairs of modes. 

However, Amblard et al. [l3| used simulations to show 
that the minimum variance quadratic estimator built out 
of the temperature field is biased if the lensing field is 
non-Gaussian at the level expected for real structure. 
Subsequently Maturi et al. @ helped to explain this fact 
by deomonstrating that the reconstructed density field 
for this estimator is biased low in regions around large 
clusters. 

In this paper, we analyze the origin of the bias and 
show that it can be nearly eliminated using the increas- 
ingly well measured properties of the CMB, and that the 
modifications we introduce produce only a very modest 
degradation in the statistical reconstruction of cluster 
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properties. Since the temperature field estimators read- 
ily generalize to polarization, these can be used as consis- 
tency checks on each other or indeed on any other cluster 
mass estimates. 

The outline of the paper is as follows. In f|TT]we review 
the general construction of quadratic estimators and the 
underlying approximations upon which they are based. 
We determine the origin of the biased estimates and 
show that they are associated with the miss-estimation of 
source gradients in the moderate to strong lensing regime. 
The bias can be eliminated at negligible cost to the signal 
to noise by imposing a strong filter against this and other 
contaminants. In ^III[ we consider the idealized perfor- 
mance of the temperature and polarization estimators. 

For illustrative purposes, we use throughout a flat 
ACDM cosmology with matter density il m — 0.24, 
baryon density fi&/i 2 — 0.0223, Hubble constant h = 0.73, 
scalar tilt n s = 0.958, amplitude eg = 0.76. 

II. QUADRATIC ESTIMATORS 

Lensing is a surface brightness conserving remapping of 
the intrinsic temperature and polarization fields from re- 
combination. Given an unlensed temperature field T(n) 
and Stokes Q(n) and £/(n) fields, where n denotes the 
angular position on the sky, the lensed fields are given by 

T(n) = T(n + V0), 
Q(n) = Q(n + V0), 
U(h) = t/(n + V<£). (1) 

Here <j> is the deflection potential, V</> is the deflection 
angle, and they are related to the convergence as 

V 2 = -2k. (2) 

All derivatives throughout are angular derivatives on the 
sky. We will furthermore employ the flat sky approxima- 
tion in this paper but all expressions readily generalize 
to the full sky with the replacement of ordinary deriva- 
tives with covariant derivatives and Fourier modes with 
spherical harmonics [13, [ID, [l6| . 

The estimators of <f> or equivalently k introduced in 
0, Q rely on two related but independent linearization 
approximations: the gradient approximation and lin- 
earization in the convergence. Only the latter, which is 
equivalent to considering lensing as a small correction to 
the source field, requires modification for use in cluster 
reconstruction. 

To see how these approximations enter, let us first con- 
sider the case of the temperature field. The analysis read- 
ily generalizes to polarization as we shall see. 

Fundamentally the temperature estimator is built out 
of a lensing induced correlation between the temperature 
field and its gradient [13] • This correlation arises from 
the gradient approximation, valid when the deflections 
are small compared with the structure in the unlensed 



field. This approximation of course does not hold for all 
Fourier modes in the CMB fields and lens (see e.g. fig). 
It need only hold the modes that are correlated by the 
reconstruction. The lensed field Tl (n) can be prefiltered 
in Fourier space to isolate modes for which the gradient 
approximation is valid 

L T (n) = T L (n), 

T L (n) = J ^Le^WfTr, (3) 

where is the Fourier filter, the subscript "L" denotes 
the lensing-filtered field and T\ is the Fourier representa- 
tion of the full temperature field. 

The lensed temperature field Tl can be approximated 
by a Taylor expansion 

T L ^f L + (Vf ■ V0) L , (4) 

which we will call the gradient approximation. For the 
case of cluster lensing, the typical deflections are < V 
compared with structure in the unlensed CMB with co- 
herence of ~ 10'. Eq. (0| is therefore an excellent ap- 
proximation for the full field even in the strong lensing 
regime. In this case the minimum variance filter can be 
used H 

W? = (C? T + ND- 1 ■ (5) 

Here Cf T is the lensed CMB power spectrum and Nf T 
is the noise power spectrum. More generally can be 
chosen to suppress modes which violate the gradient ap- 
proximation or which are contaminated by foregrounds. 

Given that the gradient approximation induces a corre- 
lation with the unlensed temperature gradient, one forms 
a quadratic estimator by multiplying Tl by a filtered gra- 
dient of the lensed temperature field 

G T (n) = VT G (n), 

Tain) = J ^L^wFTu (6) 
where Wf T is another Fourier space filter. The filter [1| 

wr=cr(c? T +Nrr i (7) 

yields a minimum variance reconstruction under the fully 
linearized approximation. Here Cf T is the unlensed 
CMB power spectrum and Eq. ([7]) is essentially a Wiener 
filter for the unlensed gradients (c.f. [6( for the more 
problematic Wiener filter for the total gradient). 

This filter has to be modified in the presence of rare 
non-Gaussian structure where lensing effects are moder- 
ate to strong. Expanding the product in the gradient 
approximation, we obtain 

G T L T » (Vf G )f L + (Vf G )(Vf • V4) L 

+[(V^)-V)Vf + (VT • V)V^)] G Ti (8) 
+[(V0) ■ V)Vf + (Vf • V)V0] G (Vf ■ V0) L . 
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The second approximation is that the quadratic term in 
(j> is negligible. If so, averaging over realizations of the 
unlensed temperature field yields a quantity proportional 
to with a proportionality coefficient related to the 
well-determined two-point function or power spectrum of 
the unlensed CMB. This quantity is therefore a quadratic 
estimator of the deflection angles V0 or k in the linear 
approximation. 

The quadratic term in <f> comes from the change in the 
temperature gradients due to lensing. Qualitatively, its 
omission is equivalent to considering lensing as a small 
perturbation to the unlensed image. Quantitatively it in- 
volves the assumption that the deflection angles are small 
compared with structures in the lens. This approxima- 
tion is violated around a cluster and the Wiener filter, 
which is based on the variance of typical regions, does 
not sufficiently suppress this region in the absence of de- 
tector noise. 



To see the effect of the quadratic term, consider the 
lensed temperature field to be filtered for small scale fluc- 
tuations as in Eq. (JSJ) with Nf T = and likewise take 
the unlensed temperature gradients VTq ~ const. In 
that case, the estimator becomes 
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FIG. 1: Gradient of the unlensed temperature T and polar- 
ization E fields as a function of the maximum I — la of a top 
hat low pass filter. For both fields, gradients are dominated 
by modes with I < 2000. 



low pass filter Wf T out to I = Iq 



t t t 



G L 



VT + (VT • V)V0) (VT • V(j)) L . (9) 



The quadratic term carries a coherent contribution whose 
strength depends on k since W 2 (j> — —2k. As k becomes 
0(1) in the interior of the cluster, the estimation will 
be biased low. In fact for k > 1 the observed gradient 
reverses and a second, flipped image of the intrinsic gra- 
dient appears at the center of an azimuthally symmetric 
cluster. 

Another way to see why the reconstruction is biased 
low in the single image regime is that the cluster mag- 
nifies the background image and decreases the observed 
temperature gradient behind it. While this bias is miti- 
gated by the Wiener filtering of Eq. ([7]), the filter cannot 
remove it entirely. Furthermore, a truly optimal filter 
would require knowledge of the cluster mass to be esti- 
mated. 

The bias arises because of the overlap in scales between 
the unlensed gradient field and the lensed temperature 
field. Maturi et al. @ suggested that it can be fur- 
ther mitigated by a combination of more strongly high 
pass filtering the lensed temperature field in Eq. ([5]) and 
utilizing the direction of the large scale gradients. The 
efficacy of this technique depends on both the assumed 
signal and the instrument noise. 

Instead let us exploit the fact that we have robust prior 
knowledge about the unlensed CMB spectrum. A hand- 
ful of well determined cosmological parameters fixes its 
shape out through the damping tail where the small scale 
gradients pick up most of their contribution. Fig.[T]shows 
the unlensed rms gradient as a function of a step function 
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Almost all of the gradient comes from I < 2000 due to 
diffusion damping of the acoustic peaks. On these scales, 
the lensing correction to the gradient is tiny for a typical 
cluster. This suggests that if we impose a sharp filter on 
the gradient to exclude higher multipoles, we lose very lit- 
tle signal and gain a clean separation between the lensed 
and unlensed structure. Such a filter is also desirable for 
large scale reconstruction to eliminate pairs of modes that 
violate the gradient approximation. Furthermore, it also 
reduces contamination from foregrounds and systematics 
that might appear as signal without a clean separation of 
scales between the two fields of the quadratic estimator. 

Now let us make these considerations concrete and gen- 
eralize them to the full set of temperature and polariza- 
tion fields X, Y £ T, E, B. E(n) and B{n) are real valued 
fields that are related to the Stokes parameters as 



Q(h)+iU(h) 



d 2 l 



[Ei + iB^e^e 11 ^ , (11) 



where ipi is the angle between 1 and the axis which de- 
fines Stokes Q. We shall assume that B is absent in the 
unlensed CMB. 

The angular space convergence estimators 



d 2 l 

(2tt) 2 



H-&-XY 



(12) 



are constructed from Fourier space estimators that are 
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quadratic in the observed fields 



\XY 



A XY 



d 2 ne- m l Re{V ■ [G XY (n)L y *(n)]} . 



(13) 

The gradient field G XY is built out of the lensed X field 
as either a spin-0 or spin-2 field depending on the spin of 
Y 



G XY = il W l 



XY X X x { 



1, 



Y = T . 



3 2i Vl y = E,B . 



L Y represents the lensed fields 



Lj = W Y Y X x { l J t . 



1, Y = T , 

' Vl , Y = E,B . 



(14) 



(15) 



For the polarization fields, they are the Stokes field Q+iU 
filtered for the E and B components. The 2 gradi- 
ent fields and 3 lensed fields produce 6 estimators from 
the temperature and polarization fields for consistency 
checks. In addition, replacement of the divergence in 
Eq. (fT3"|) with curl should leave estimators that are con- 
sistent with noise. 

A x Y is a normalization coefficient set to return unbi- 
ased estimators under the fully linearized approximation. 
For arbitrary filter functions, it is determined by noting 
that @ 

(X h Y h )=2f XY (h,l 2 )K l , (16) 
where 1 = L + I2 ^ and 

f TT (h,h) = (7^(1.10 + ^(1-12), 
f TE (h,h) = C? E cos 2 Vllla (l- h) + Cl E (l-h), 
f TB (h,h) = C£ E S m2<p lll2 (hh), 
f ET (h,h) = f TE (h,h), 

f EE (h,h) = [C EE (l-h) + C EE (l-h)}co S 2^ li2 , 
f EB (h,h) = C EE sm2<p llh (Yh). (17) 

Here y'lib = fii — fla- ^he normalization is given in 
terms of these quantities and the filters as 

where 

c T (li,l 2 ) = 1, 

C B (li,l 2 ) = 008 2(^1^2 , 

c B (h,h) = sm2p llh . (18) 
For the lensed field filter, we retain the choice of [1,0] 



w y = (c, 



YY 



N, 



YY\-1 



(19) 



As described above, the first modification is that the gra- 
dient weights are set to zero above Iq = 2000 



Fig. Q] shows that this is also a good choice for the E field. 
For I < Iq, we retain the Wiener filter 



W, 



XY _ qXY triXX 



N, 



XX\-1 



(1<Ig), (21) 



W, 



XY 



0. 



(I > Ig) 



(20) 



for Y ET,E and 

w xb = c XE {C xx + iV^^)" 1 , (/ < l G ) , (22) 
for Y = B. 

Secondly, we distinguish between TE and ET estima- 
tors. With the gradient filter employing only I < Iq, 
the scale symmetry between the two fields is broken. It 
is then advantageous to separate the estimators. The 
TE estimator uses T for the gradient field and E for 
the lensed field. If contamination from unpolarized clus- 
ter emission is strong then this estimator can be used to 
eliminate its effects. The ET estimator uses E for the 
gradient field and T for the lensed field. This estimator 
places less demanding requirements on the experiments 
in that only the relatively large E field of the acoustic 
peaks need be measured to high signal to noise. Indeed 
this measurement need not come from the same experi- 
ment that measures the cluster signal. Both estimators 
strongly reject contaminants from foregrounds and sys- 
tematics. Only the background CMB will have T and 
E correlated and anticorrelated in the specific oscillatory 
pattern of the acoustic peaks. Unfortunately, the sam- 
ple variance of these estimators is also higher given the 
imperfect correlation between the two fields. 

On the other hand, the sample variance of the EB es- 
timator is reduced by the fact that B modes are assumed 
absent in the unlensed fields |T(3] • In practice the noise 
of the actual estimator will depend on the B-modes con- 
tributed by foregrounds and other contaminants. Like 
the ET estimator, cross correlation of the large scale E 
gradients with the small scale _B-modes should help re- 
duce any bias arising from contaminants. Since the ex- 
perimental requirements for EB are similar to that of 
TE, TB, and EE but yield better prospects for con- 
straints, we focus on the EB, ET and TT estimators in 
the next section. 



III. IDEALIZED EXAMPLES 

To illustrate the performance of the estimators, let us 
take the idealization that all lenses are NFW [19( dark 
matter halos and the observed CMB fields have no con- 
taminants aside from white detector noise. We will ad- 
dress more realistic cases in a future work [20]. The 
estimators should remain nearly unbiased since to pass 
through the reconstruction filters, a contaminant must 
be antisymmetric around the cluster center with an axis 
whose direction is correlated with the large scale gradi- 
ents. In addition, contaminants will not appear in the 6 
quadratic estimators, especially those involving cross cor- 
relation, in the same way. They can however contribute 
substantially to the noise and we will examine here the 
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performance of the estimators as a function of detector 
noise as a proxy for other contaminants. 
The NFW density profile is given by [n| 



p(r) 



1 



r/r s (l + r/r s ) 2 



(23) 



where the normalization coefficient can be expressed in 
terms of the mass of the halo. We define the mass to 
be that enclosed at riso defined to be the radius which 
encloses the mass at an overdensity of 180 times the mean 
density 



riso 



3M 



1 



47T 180tt m p c 



1/3 



(24) 



where p c is the critical density today. Likewise, we define 
the concentration of the halo in terms of this radius 



ciao — 



^180 



The convergence profile for such a halo is [2l| 

= -t-HqD l ——{1 + z L ) T — 

4tt D s p c rj /(ciso) 



(25) 



(26) 



where the projected scale radius 9 S = r s j Dl. D is the 
comoving distance in a flat universe with subscripts "L" 
denoting the distance to the lens, "S" denoting distance 
to the source, and "LS" the distance between the lens 
and source. A primary advantage of the CMB is that 
the source is behind all clusters and at a well-determined 
distance Ds — 14.3Gpc. The concentration factor is 



/(c) = ln(l + c) 



1 + c 



and the functional form of the profile is [21] 



9{x) 



1 



: atan 



1 



x 2 - 1 
1 
3 



1 



2 , l-x 

zatanhi / 

1+x 



(X = l). 

For illustrative purposes, we take 



Mi 



M 



= 2. 



(27) 

(x>l) 

, (x < 1) 
(28) 

(29) 



ciso = 3.2 and Zl = 0.7. For these parameters 8 S = 0.94' 
and k(8 s ) = 0.1. Note the low value of the concentration 
reflects the low as fiducial cosmology and is determined 
by the scaling of [22|]. In a higher erg cosmology, clusters 
at a fixed mass are more concentrated but correspond- 
ingly larger clusters are more abundant. Other numbers 
are typical of SZ selected clusters with upcoming surveys. 




FIG. 2: Temperature (TT) reconstruction in the absence of 
detector noise and other contaminants. From left to right, top 
to bottom the reconstructions are stacked with 1, 10, 10 2 , 10 3 



identical clusters of Mu = M/lO 14 ^ _1 Af = 2 and z 
Shown here is a 35' x 35' region around the cluster. 



0.7. 



For the simulated reconstruction we consider a fiducial 
experiment with a #fwhm = 1' beam and varying levels 
of noise on the sky 



N, 



TT 



= AiB 



N t BB = A 2 P B- 2 . 



where the beam factor is 
Bi = e" 



2 fFWHM/161n2 



(30) 



(31) 



with #fwhm in radians. It exponentiates the white de- 
tector noise on the beam deconvolved sky. We further 
assume that Ap = ^/2At- For this fiducial beam, a 
pixel scale of 9 P i X = 0.2' suffices. The field is chosen to 
be 512 x 512 pixels (100' x 100'). 

We lens realizations of Cf T , Gf E according to the full 
non-linear prescription of Eq. |1| with 6th order polyno- 
mial interpolation of the source image. This procedure 
allows strong lensing effects to appear at the center of 
the cluster. To this lensed image we add a realization of 
the noise power spectra above. 

In Fig. [21 we show the reconstruction with At = sim- 
ulations of the TT-estimator for N = 1, 10, 10 2 , and 10 3 
stacked clusters. In this detector noise free limit, there 
is a clear detection with only one cluster. It appears as 
a concentrated spike on top of a noisy but slowly vary- 
ing background from chance correlations of the Gaussian 
random background field. The reconstructed single im- 
age lacks the azimuthal symmetry of the lens but symme- 
try is recovered by stacking images 0]. This asymmetry 
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comes from the fact that the deflections cannot be esti- 
mated in the direction orthogonal to the gradient. The 
estimator however is unbiased once many realizations of 
the gradient are stacked. 
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FIG. 3: Recovered cluster-mass correlation function or 
stacked profile with 24000 detector noise free simulations. 
Bands represent the rms fluctuation per fOOO clusters in the 
sample. For M\a = 2, the reconstruction has no detectable 
bias at the I — 2% level within a few arcmin. 

Stacking is equivalent to measuring the cluster- 
convergence correlation function 

U(0) = /^<*(M)>, (32) 

where the coordinates are centered at the location of the 
cluster such that 9 is the separation from the center and <f> 
is the azimuthal angle around the cluster. The averaging 
is over the clusters in the sample. By defining the observ- 
able as the cluster-convergence correlation function, the 
estimator naturally generalizes to non-identical clusters 
and projected mass along the line of sight associated with 
the cluster. In practice, we evaluate the correlation func- 
tion at discrete intervals in pixel units by interpolation 
on the stacked image. Neighboring bins out to several 
arcminutes are therefore highly correlated as can be seen 
directly in Fig. [2j This correlation is in fact useful for 
distinguishing the signal from noise. 

In Fig. [31 we show the recovered correlation function 
from 24000 detector noise free realizations along with 
the scatter per 1000 clusters with M\& = 2. We have 
smoothed both the lens and reconstruction with a 1.5 
pixel FWHM beam after the fact for comparison. The 
estimator has no detectable bias at the ~ 1 — 2% level 
in the well-measured regime within a few scale radii. At 
five times the fiducial mass or Mi 4 = 10, a low bias 
of ~ 8 — 9% develops for the fiducial gradient cut of 
£<3 = 2000. For these very rare clusters, the signal is 
large enough that a more aggressive Iq = 1000 yields 
a bias free reconstruction with sufficient signal to noise. 
Alternately, the bias can be calibrated out with simula- 
tions. 
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FIG. 4: Beam filtered correlation function from 12000 clus- 
ters. The reconstruction is low pass filtered at I < £beam(= 
8095) to remove contributions below the beam scale. Shown 
for comparison are the filtered input lens of M14 = 2 and a 
trial model of M14 = 2.6 for comparison. With a noise level 
of 10/iK' the two cases are clearly distinguished. 

For cases with finite noise, we explicitly remove effects 
below the beam scale for numerical convenience. We low 
pass filter both the lens a nd th e reconstruction with a 
step function at Zbearo — V8 In 2/ #fwhm ■ This also has 
the desirable effect of making the reconstruction less sen- 
sitive to potential contamination near the cluster center. 
The filtered correlation function is shown in Fig. [4] and re- 
mains unbiased and well-defined. It is equivalent to mea- 
suring the cross power spectrum at multipoles I < Zbeam- 
Fig. [4] also shows that with A T = 10/xK', Mu = 2 and 
2.6 are significantly distinguished. 

To quantify this sensitivity, we compute the covariance 
matrix of the correlation function from the simulations 

0(6,6') = (ie K (6)U(0')) -UWUW) > (33) 

where the discretization again is in bins of the pixel 
width. To estimate the significance with which two mod- 
els can be separated, we evaluate A% 2 between a test 
model £,' CK {0) and the true model £ C k(0) 

ax 2 = - umc-^cv) - e CK (0')] , (34) 

e,G' 

summed out to 6 = 5'. Beyond a few scale radii of 
the cluster the correlation function is dominated by large 
scale structure. 

In Tab. U we show the detection significance, the dif- 
ference between the fiducial M14 = 2 lens and zero with 
the covariance matrix evaluated at zero signal. We have 
scaled the significance per cluster by VN for N = 10 3 . 
At a noise level of 10/jK', the significance or signal to 
noise is y ^X 2 ~ 0-4 per cluster. Above this noise level, 
only TT and ET can provide significant detections. Al- 
though TT has roughly three times the signal-to-noise, 



7 




1 10 



Aj=Ap/2 1/2 (nK') 

FIG. 5: Idealized mass sensitivity for #fwhm = 1', 9 < 5'. 
The fractional errors on the idealized NFW mass with 1000 
clusters for the TT, ET and EB estimators. Unrelated mass 
along the line of sight will add noise to the estimates as will 
any emission from the cluster. No information from the re- 
construction is used below the beam scale I > /beam- 

ET can still provide a useful check against contamination 
and systematics for a first detection. Recall that the filter 
construction involves the oscillatory acoustic TE correla- 
tion and strongly rejects spurious signals. At noise levels 
between 1 — 3/iK', the EB estimator has both the signal- 
to-noise and the ability to reject contaminants to make 
it competitive with TT. Below 1/tK', the EB estimator 
dominates. 

In Fig. [5] we plot the sensitivity to the cluster mass. 
Specifically we calculate the fractional change in mass 
that would generate a A\ 2 — 1 with 10 3 clusters. The 
covariance is evaluated at the true model Mn — 2 so as 
to include the sample variance of the unlensed gradients. 
The same features found for detection significance hold 
for mass sensitivity. For example at 10/xK', TT, ET, and 
EB provide 13%, 43%, 53% mass estimates whereas at 
lfjK' the numbers improve to 5%, 14%, 6%. 

Finally, these sensitivities depend on the beam scale 
only weakly near the fiducial #fwhm = 1' and = 
1 — 10/iK' level. For example, at 3/xK', improving the 
beam to 0.5' has negligible effect on the signal to noise. 
Enlarging the beam to 2', degrades the TT and ET esti- 
mators by ~ 30 — 40% in mass sensitivity. Since the EB 
estimator is noise dominated at the 1' scale even without 
beam, it degrades negligibly. 

TABLE I: Detection significance y/Ax 2 or S/N for 10 3 clus- 
ters of M14 = 2, z = 0.7 and #fwhm = 1'. 



Type 


3/iK' 10/xK' 30/tK' 


TT 


22.9 11.5 6.1 


ET 


8.0 3.5 1.4 


EB 


8.6 2.4 0.6 



IV. DISCUSSION 

We have developed mass estimators from 6 quadratic 
combinations of CMB temperature and polarization 
fields that retain their minimum variance characteristic 
for large scale structure and provide nearly unbiased re- 
constructions for rare clusters where lensing effects are 
moderate to strong. The key difference from previous 
work [8, 9] is that the Wiener filtering for unlensed CMB 
gradients is augmented with a sharp filter that removes 
any gradients at Iq > 2000 that are not part of the source 
fields from recombination. 

This sharp filter should also help prevent false signals 
from cluster emission such as the thermal and kinetic 
Sunyaev-Zel'dovich effect, point source emission as well 
as other foregrounds. To even appear as excess noise, the 
contaminant at ~ 1' must have a component that is anti- 
symmetric about the cluster center. To bias the measure- 
ment, the direction of the asymmetry must be correlated 
with both the polarized and unpolarized contamination 
at ~ 10' in the same way as the acoustic oscillations in the 
CMB. Note that projection effects from mass associated 
with the cluster should be considered part of the signal 
and can be calibrated by iV-body simulations. True con- 
taminants however can substantially increase the noise of 
the estimators. An evaluation of the efficacy of the filter 
in the presence of real world contaminants is beyond the 
scop e of this paper but will be treated in future work 

m 

Here we have tested the estimators against idealized 
signal and noise: identical clusters lenses with NFW dark 
matter profiles in the presence of beam filtered white 
noise. At a noise level of 10/iK' and with a 1' beam, 
the TT estimator can provide a ~ 10cr detection per 10 3 
clusters of M ~ 2 x 1O 14 /j _1 M or equivalently a ~ 10% 
average mass measurement. The ET estimator, based on 
a separate measurement of _E-polarization in the acous- 
tic peaks and T measurements around the cluster, can 
provide an important cross check on a first detection. 
Since the estimator filters for the TE correlation and an- 
ticorrelation of the acoustic peaks it provides a strong 
discriminant against contamination and systematics at 
the price of a factor of ~ 3 in signal-to-noise. Ultimately 
with experiments that produce foreground free maps of 
the polarization to < 3/iK', the EB estimator will re- 
turn the best estimates, with the potential to ultimately 
provide ~ 1% measurements. 
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